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Abstract 

The mass function of primordial black holes created through the near-critical gravitational collapse 
is calculated in a manner fairly independent of the statistical distribution of underlying density 
fluctuation, assuming that it has a sharp peak on a specific scale. Comparing it with various 
cosmological constraints on their mass spectrum, some newly excluded range is found in the 
volume fraction of the region collapsing into black holes as a function of the horizon mass. 
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In the study of gravitational collapse with formation of a black hole a critical phenomenon 
was discovered by Choptuik [|J , who performed a series of numerical calculation of the evolution 
of a spherically symmetric system with a minimally coupled scalar field. He showed that the 
resultant mass of the black hole, if formed at all, scaled as Mbh oc (p — p c ) 7 with p > p c , where p 
is a control parameter of the system and p c is the critical value above which a black hole forms. 
The remarkable feature of the above formula is that the power-law exponent was found to be 
independent of the choice of the control parameter. Since then a number of authors investigated 
various systems with scale-free matter ingredients and found essentially the same result with the 
unique value of the power index, 7 ~ 0.36 HQ. An analytic explanation has also been given by 
Koike et al. Although this is certainly an interesting subject of study in general relativity, its 
relevance to the black holes created in astrophysical or cosmological setting had not been known. 
Recently, however, Niemeyer and Jedamzik || observed the same type of the critical behavior 
in their new numerical calculation of the formation of primordial black holes (PBHs) out of 
radiation fluid in the Friedmann-Robertson- Walker background. 

PBHs are formed in the radiation dominated era when a perturbed region enters the Hubble 
radius if the amplitude of the fluctuation exceeds some critical value (5 C ~ 1/3 ||, which has been 
confirmed by numerical analysis done two decades ago pyTcfl. It has also been concluded that 
the resultant PBH has a mass of order of the horizon mass at formation. Consequently all the 
studies done so far to constrain the mass spectrum of PBHs are based on the assumption that 
PBHs of a specific mass is created at a specific epoch, namely, when the mass scale entered the 
Hubble radius as long as we are concerned with those formed during radiation domination [11]. 



The same assumption has also been made in building models of inflation to realize formation of 
PBHs on some specific mass scales p| 111. 



According to Niemeyer and Jedamzik ||, however, PBHs are not only produced with about 
the horizon mass but also with much smaller masses with the scaling formula 

M BH = K(5 - 5 C )\ (1) 

where as the control parameter 5 they chose the additional mass in the perturbed region in unit 
of the horizon mass when the relevant scale entered the Hubble radius. They investigated two 
classes of configurations of the perturbed region and found a universal power index of 7 = 0.35 
in agreement with Evans and Coleman || and Koike etal. Q. They have also confirmed it is 
independent of the choice of the control parameter || , as it should be. 

Although the scaling relation (|l|) is expected to be valid only in the immediate neighborhood of 
S c , most black holes are expected to form with such an initial value of 5, because it has generically 
a rapidly declining probability distribution function (PDF) near 5 = 8 C and the probability to find 
5 3> S c is exponentially smaller. Hence it is sensible to calculate the expected mass function of 
PBHs using the formula (|l]). The above-mentioned feature of the PDF also allows us to estimate 
the mass function fairly independent of the specific form of the PDF of primordial density or 
curvature fluctuations as seen below. 

The purpose of the present paper is to investigate the effect of the near-critical gravitational 
collapse on the cosmological constraints on the mass spectrum of PBHs and on the model building 
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of inflation to produce PBHs. As a result we find that, thanks to the fact that 7 is relatively 
small, we do not have to change the conventional view to the issue drastically except that we find 
a new constraint on the fraction of the space collapsing to PBHs when the horizon mass is in the 
range between 4 x 10 14 g and 6 x 10 16 g. 

In this paper we study the case of a primordial spectrum of density fluctuation that is sharply 
peaked on a single but arbitrary mass scale just as in the initial configuration adopted in 1|. 
Several models have been proposed to realize such a spectral shape in inflationary cosmology 
|l2,15-18|]. Even in the case with such a simple spectrum, the configuration of a perturbed region 



has functional degrees of freedom and one must in principle calculate both the probability to 
realize each initial configuration and the resultant mass of the black hole, if formed at all, in 
order to calculate the mass function of PBHs. In the language of the near-critical collapse, there 
are many possible one-parameter families of the initial data in the configuration functional space 
to approach the critical surface and we find different values of 5 C and K for each family with the 
exponent 7 being the only universal quantity. Hence we should first calculate the mass function of 
PBHs when the initial configuration is changed along each one-parameter family near the critical 
surface and then integrate over such families in the functional space together with their relative 
probability to obtain the overall mass function, which is a formidable task. Fortunately, however, 
there is yet another common feature in PBH formation, which is the fact that the typical mass 
of the black hole is of order of the horizon mass independent of the shape of the perturbed region 
even in the presence of the critical phenomenon as shown explicitly in ||. We thus reduce the 
problem with infinitely many degrees of freedom to a one-dimensional issue, that is, we adopt an 
approximation that the system is governed by a single parameter, 5, characterizing the amplitude 
of fluctuations, and that the most common black hole has the horizon mass, Mh, when the peak 
of fluctuation entered the Hubble radius at t = tn- 

Since the PDF of 5, P(S), is a steeply declining function around 5 C we write it as 

P(8)dS = e~ m d5, (2) 

where f(5) is a well-behaved function around 5 ~ S c . If 5 is Gaussian distributed, as assumed in 
most literature fj,19,20fl, f(S) is explicitly given by 



f(5) = ln(v&) + (3) 

with a being the dispersion of 5. 

Then the probability, (3(Mh), that the relevant mass scale has an above-threshold amplitude 
of fluctuations to collapse into a black hole as it enters the Hubble radius is given by 

f3(M H ) = I P(S)d5 = f e- f Wdd, (4) 

which is also equal to the volume fraction of the region collapsing to a black hole at t = tu ■ Due 
to the fact that P{5) is a steeply decreasing function the integral is sensitive to only its lower 
bound 5 C . Furthermore we can Taylor expand f(5) as 
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to find 



f(S) = f(5 c ) + f'(S c )(5 -S c ) + --- = f c + q(S-6 c ) + ---, 



P(M H ) 



(5) 
(6) 



The above formula (g) based on the approximation (|B|) is applicable if \f"(S c )\ <g 2 . In the case 
of a Gaussian distribution this corresponds to a <C <5 C . Putting <5 C = 1/3 and using @ in @ 
we recover Carr's formula ||], f3 ~ <r exp[— 1/(18<7 2 )]. If P(5) is non-Gaussian, we must analyze 
case by case, but as long as we consider the case with small enough (3{Mh) in the cosmologically 
allowed range jnj and use specific non-Gaussian distributions reported in the literature fi"6|,18,2Tf , 
we may justify (||). Hence this simple approximation enjoys a wide applicability. 

Then the contribution of PBHs to the density parameter at formation 1 1 zltn, is given by 



1 



with S(M) 
reads 



M H 

5 C + (M/K)^ 



M BU (5)P(5)d5 



1 



(7) 



Using (H) and 
dQ BU (M, t f ) 



M H jo 

in the above integral, the differential mass spectrum 



dlnM 
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~K 
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M 
~K 
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which has a peak at M max = Kq 
stated above, let us identify M ma 
then find 



7 (1 + 7) 7 . Since the typical black hole mass is around Mh as 
K with Mh following the line of thought explained above. We 
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which is now independent of q and K. Formal integration of @ from M = to oo yields 
Q BU (tf) ~ (1 + 7)~ 1 r(7)/3(M H ) = 0.80/9(Mh) for 7 = 0.35, as opposed to the case all the PBHs 
have the horizon mass where we would find ^bh(^/) = /3(Mjj). Thus the total mass density 
of PBHs does not change practically even if we take the effect of the near-critical collapse into 
account. 

The mass spectrum (||) is depicted in Fig. 1 where we find that the abundance of smaller- 
mass black holes are suppressed even in the presence of the critical behavior. We can hence 
convince ourselves that the previous assumption that only horizon-mass black holes are likely be 
produced is a good approximation in model building where we try to identify some astrophysical 
objects such as MACHOs [§2| with PBHs H©®. Indeed ' if to ° man y Jupiter-mass PBHs 
were produced simultaneously in an attempt to produce MACHO-mass PBHs, PBH-explanation 
of MACHOs would be ruled out because the observation of the MACHO group is very sensitive 
to Jupiter-mass objects and their abundance in the galactic halo has already been stringently 
constrained |23|. However, the mass fraction (|9|) implies we are free from such a problem. 
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FIG. 1. Expected mass function of PBHs (|9|) created through the near-critical gravitational 
collapse in the case density fluctuation has a sharp peak on the mass scale Mh- 

We must, on the other hand, consider its effect on cosmological constraints more carefully. At 
later time t we find the mass spectrum, !F(M, t;tf), of the holes produced at tf as 



dn BR {M,t;t f ) = cin BR (M, t f ) ft 
dlnM dlnM I 



d^BH(M, t f ) ft 
dlnM ' 




where 



t eq = 4.2 x 10 iU (f} /i ) sec = 1.9 x lO^r^sec, T = 



Q / h 
03 V 07 



2 On/l 2 



0.147' 



(10) 



(11) 



is the equality time with Qq and h being current values of the total density parameter and the 
Hubble parameter in unit of lOOkm/sec/Mpc, respectively. 

Cosmological constraints on the mass spectrum of the PBHs can be classified to three classes. 
The first one applies to heavy black holes with mass Mbh > 4 x 10 14 g which have not evaporated 
by now. Current mass density of such holes should not exceed the total mass density of the 
universe. The second class is due to the radiation of high energy particles from evaporating 
black holes p4] , p5| . Various constraints have been imposed from primordial nucleosynthesis |26j, 
microwave background radiation [27|, and gamma-ray background radiation pq . Finally, there 
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may be yet another class of constraints if PBHs do not evaporate completely but leave relics with 
mass of order of the Planck mass or larger |2{|. Their mass density should remain small enough. 

Even if we adopt the new mass spectrum ([]]), cosmological constraints of the first and the 
third classes described above are not altered practically because the total number of black holes 
as well as the abundance of most typical black holes, which have approximately the horizon mass 
at formation and dominates their mass density, are not expected to change significantly as shown 
above. In particular, the constraint of the first class, namely, that from the total mass density of 
PBHs reads 



P(M H ) < 2 x l(T 19 r 



M H 



4 x 10 14 g 



0.3 : 



(12) 



demanding that Qbh today should not exceed Qq. 

We now concentrate on the constraints associated with evaporation. The lifetime of a black 
hole with mass M is given by 

/ M \ 3 

t ev (M) = 1.7 x lO 3 ^ 1 (j^-J tpi, (13) 
where Mpi (tpz) is the Planck mass (time) and g is the effective number of massless state radiated 



from the black hole in units of 7.25 which applies at the low-energy or large mass limit [20]. 
A black hole will evaporate before the equality time if its mass satisfies the inequality M < 
6 x 10 12 <j/3r~3g = M eq . Eliminating tf from both sides of ( |l0|) by virtue of the relation Mp = 
M 2 Pl t H ~ M Pl tf, we find 
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The above mass spectrum should be compared with the following constraints of the second 
class [0JT|j2|,|25|-|8|. 



F{M,t ev (M)) < 1 x 10- 2 M^ , 

< 5 x 10- 7 m!, 

< 2 x 1Q- 8 M? , 



< 
< 



4 x 10" 2 , 
1 x 10" 8 , 
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for M 
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6 x 10 n g, 



(15) 



6 x 10 



ii 



10 13 g 



for M ~ 4 x 10 14 g, 



ao 



where M 10 = M/10 1 

Putting M = Mp in (|14|) and comparing it with (|T5|), we practically recover the constraints 
on (3(Mp) previously obtained in the literature, which are depicted by a solid line in Fig. 2. In 
the presence of the near-critical collapse, however, we should also consider constraints on (3{Mp) 
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arising from PBHs with mass M <C Mh- In fact, however, due to the steep shape of ( |l4| ) we 
find a newly excluded range Q only in the region Mjj > 4 x 10 14 g, which is imposed from the last 
constraint of (|l5|), namely, 
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FIG. 2. Cosmological constraints on the volume fraction of the region, (3, above the threshold 
of PBH formation when the horizon mass is equal to Mh- The region between the solid and the 
dotted lines is newly excluded. 

As is seen in Fig. 2, this new constraint is more stringent than ( |l~2| ) for Mh < 6 x 10 16 g. 

It has been argued that high energy phenomenon associated with currently evaporating black 
holes with mass ~ 4 x 10 14 g can explain the origin of a class of gamma-ray burst [pCfl . If this is 
the case, their abundance should be around f^BH = 10 -8 today. If such a tiny amount of black 
holes were created at the low-mass tail of the near-critical collapse, one could explain the origin 
of both such bursts and dark matter simultaneously. That is, PBHs with mass around 6 x 10 16 g 
could be the dominant dark matter component which makes up 0, ~ 0.3 today. 



There may be another small newly excluded region around Mh ^ 10 g. However, we did not 
depict it in Fig. 2 since it is based on a rather qualitative constraint from nucleosynthesis obtained 
two decades ago, and refined calculations may well alter it. 
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In summary, we have obtained an approximate but generic form of the mass function of 
primordial black holes which are produced through near-critical gravitational collapse of radiation 
fluid in the early universe in the case density fluctuations have a sharp peak on a specific scale. It 
is fairly independent of the statistical distribution of density fluctuations. Due to the smallness of 
the critical exponent 7 = 0.35 the resultant mass function has a steep spectrum in the low-mass 
tail, so that the previous assumption that the PBHs are created with nearly the horizon mass 
at formation is basically correct. Nevertheless we find some newly excluded range in the volume 
fraction of the region with 5 > 8 C as a function of the horizon mass as depicted in Fig. 2. 

The author is grateful to Professor Andrei Linde for his hospitality at Stanford University, 
where this work was done. This work was partially supported by the Monbusho. 
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